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东亚 飞 蝗 中 长 期 数量 预测 的 研究 
Bue Jen FRB | 


《中 国 科学 院 动物 研究 所 7 


摘要 本 文 以 洪 泽 湖 蝗 区 为 例 ， 分 析 了 东亚 飞 蝗 中 长 期 数 基 预 测 方程 的 建立 基础 ， 提 出 三 种 预测 方 
法 《一 ) 根 据 种 全 动态 型 趋势 进行 外 推 估 值 ,( 二 ) 应 用 随机 序列 及 周期 方程 的 预测 法 ,《 三 ) 多 因素 过 谊 台 妇 
预测 法 ,此 三 种 方法 虽 各 有 其 独立 的 预测 模式 ,但 在 应 用 时 应 将 其 视 为 总 预测 式 中 的 三 个 相互 补充 的 部 分 。 


”昆虫 的 发 生 预 测 工 作 , 可 以 说 自 有 效 温 积 法 则 提出 后 不 久 就 已 开始 ,但 只 限于 发 生 期 

及 发 生 世 代 方 面 ,数量 癫 测 的 研究 比较 述 ,至 于 应 用 数理 统计 方法 把 若干 主要 的 作用 因 来 
包括 在 一 个 预测 式 或 种 到 数学 模式 中 , 则 是 近 几 年 方 开始 的 工作 。 但 发 展 较 快 ,和 渤 今 见于 
文献 中 WEAH Nomadacris septemfasciata) (Symmons, 1959) ARH “sR RGA 
次 , 1959), KI SBE Diprion hercyniae), 42765 (Choristoneura fumiferana) (Mortis, 
1963), ÆR (Plutella maculipennis Curtis), RAH ERC Archips argyrospila) (Watt, 1963) 
BERKIRCOstrinia nubilalis) (Barlow, 1963), NAJE EAE (Sea 
归 )、 微 分 方程 (Logistic HAA BRIER HAA) Ri AES 

在 飞 蝗 发 生 数量 预测 方面 ， 不 少 地 区 开始 都 系 直 搂 以 去 秋 查 卵 或 查 成 虫 的 结果 为 依 
据 , 佑 计 来 年 春季 发 生 数量 (美国 .加拿大 ,1938; 苏 联 ,1950; 西 班 牙 ,1957)，, 继 而 根据 前 一 
NK AAR, AM Cet CHA RASMIER, RAS 1946 年 起 采用 此 法 预测 非洲 
地 区 的 沙漠 蝗 (Gunn，1960)， 苏 联 鲁 根据 浴 水 泛滥 高 度 预测 亚洲 飞 蝗 的 发 生 面积 〈《Caxa- 
paB，IHyMakoB8，1950), 浊 年 在 赞比亚 则 进展 到 根据 湖水 位 及 降水 量 所 计算 的 角 归 式 预测 
红星 下 一 季 的 发 生 级 数 (Symmons，1959), 而 在 澳大利亚 条 以 上 代 发 生 级 数 与 前 季 的 降水 
FETE SS ~~ RAKEI BB AY Be AE Be (Casimir, 1962). 

我 国 对 于 东亚 飞 煌 的 预测 工作 开始 于 1952 E, ROI RAKE KB RASS 
卵 的 结果 估计 来 年 春季 的 发 生 面积 , 1958 及 1959 年 先后 以 降水 景 ,温度 ,晴雨 日 数 ,湖水 
位 及 虫口 基数 为 依据 , 提出 发 生 程度 的 经 验 公式 CEER, 1958) KIEK RÆ MR 
ERBACH, 1960), 

本 文 所 报导 的 系 继续 上 项 研究 所 获得 的 近期 结果 ， 落 重 在 数值 预测 方面 。 在 分 析 方 
法 上 得 到 本 院 计 算 技 术 研 究 所 、 中 央 气 人 象 研究 所 长 测 组 ,北京 大 学 数学 力学 季 概 率 论 教研 
组 ,本 院 地 球 物理 研究 所 杨 鉴 初 教授 的 热情 帮助 ,计算 所 代为 计算 了 大 量 的 资料 ， 本 次 初 
稿 又 承 计算 所 徐 钟 济 教授 审阅 产 对 文中 数理 统计 方面 的 处 理 提出 宝贵 的 意见 ， 此 外 本 空 
陈 玉 平 同志 贿 加 了 部 分 计算 工作 , 刘 萌 增 同 志 代 为 绘图 ， 泗 洪 蝗虫 防治 站 供给 部 分 资料 ， 
一 并 在 此 志 痪 。 


1) ACHE 1960 年 昆虫 学 会 年 会 上 宣读 。 
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材料 整理 ”我 国 飞 拍 蝗 区 计 有 四 个 类 型 ,影响 飞 蝗 大 发 生 的 因素 随 蝗 区 类 型 而 不 同 ， 
本 广 仅 以 洪 泽 潮 蝗 区 为 例 , 希 望 通过 对 此 类 典型 蝗 区 飞 蝗 上 发生 因素 的 分 析 , 以 及 建立 预测 
式 的 过 程 说 明 对 昆虫 中 长 期 数量 预测 的 观点 ,方法 及 其 理论 依据 。 

纵 泽 潮 昕 区 系统 资料 的 积累 开始 于 解放 后 的 第 三 年 ,到 1963 年 止 只 有 十 二 年 资料 ， 
时 间 短 , 历史 特征 表现 不 出 来 , 因而 在 分 析 历 史 规 律 时 引用 解放 前 分 散 于 各 方面 的 资料 ， 
其 中 除 蝗虫 资料 外 ,并 以 淮河 流域 下 游 1663 一 1962 年 的 早 涝 资料 作为 分 析 时 参考 。 以 本 
区 有 系统 气象 资料 及 翔 水 位 记载 的 近 50 年 (1913 一 1962) 的 飞 凰 发 生 资料 作为 分 析 依 据 ， 
FEVA.1952—1963 年 资料 作为 检验 的 依据 。 

确定 飞 峰 发 生 指 标 时 ,考虑 到 解放 前 资料 仅 有 猩 兰 范围 轻重 程度 及 粗略 的 虫口 发 生 
情况 的 记载 ,无 具体 的 面积 及 密度 数值 ,解放 后 经 过 几 年 的 项 剂 防治 与 政 造 蝗 区 后 ， 大 面 
积 发 生 密 度 已 大 幅度 下 降 , 仅 在 未 经 防治 或 渤 后 的 局 部 小 面积 内 有 高 密度 蝗 到 的 出 现 , 若 
单独 以 平均 密度 作为 区 分 发 生 程 度 的 指标 ， 则 不 能 代表 实际 情况 ,因此 ,条 用 了 发 生 面积 ， 
与 密度 相 结 合 的 分 级 方法 ， 后 者 又 包括 大 面积 的 平均 密度 及 平均 最 高 密度 。 计 算 级 别 所 
用 公式 
了 一 互 |3 一 5 
Op Os 


虫口 密度 / 米 2 S-— REAP, onos 


D = 








式 中 DD 一 一 级 数 (面积 与 密度 的 总 和 ), P 
两 者 的 标准 离 盖 ,五 .一 二 者 的 均 数 。 
用 此 公式 分 别 计算 各 年 夏 、 秋 蝗 级 数 ( 夏 、 秋 蝗 发 生 级 数列 为 1、.2、3、4、5 级) 及 年 
发 生 级 数 ( 列 为 10 级 )。 
作用 因素 的 分 析 及 筛选 
-作用 湖区 飞 蝗 数 量变 动 的 因素 有 三 方面 ， 1) 自然 地 理 因素 , 主 要 是 气候 及 潮水 位 变 
化 ; 2) 人 为 因素 ,包括 菇 剂 防治 及 耕作 等 人 类 活动 ; 3) 飞 蝗 种 考 自 身 的 若干 生物 学 因 束 。 
根据 以 往 的 分 析 ( 马 世 骏 ,1958), 洪 举 湖 星 区 及 淮河 下 游 地 区 飞 蝗 大 发 生 的 间隔 年 
数 , 有 2 一 3 年 、9 一 11 年 不 等 , 无 明显 的 周期 性 规律 , 但 各 次 大 发 生年 距 之 间 是 否 有 一 定 
的 关联 性 ， 抑 或 属于 各 年 辣 独立 变量 的 出 现 ? 系 探索 飞 蝗 历史 独 儿 特 性 首先 应 明 稍 的 问 
题 ,作者 取 表 1 中 1913—1962 年 的 资料 ,利用 游程 公式 | 
E= —A— 20g | | 
Vapa — 3pg) a 
进行 检验 , n = 50, 一 游程 总 数 ，p 一 AREE MUUE PI a 所 得 
结果 如 下 
28 一 100(0.2016) 
EO 2016)(1 — 3 X 0.2016) 
{E p = 0.02 水 平 上 显著 。 说 明 各 年 间 的 发 生 有 一 定 规律 性 的 关联 。 
什么 是 影响 此 种 规律 的 因素 ， 以 往 分 析 表 明 飞 蝗 猛 狐 与 当地 区 的 星 浇 发 生 有 直接 关 


= 1.964 


Ro 
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Hl 洪 泽 湖 昕 区 1663—1962 年 300 FK BM RETR 
水 年 时 年 蝗虫 大 发 生年 
年 别 | 年 别 | 年 别 | 年 别 j 年 别 | 年 别 | 年 别 | 年 别 | 年 别 | 年 别 | 年 别 


1665 1736 1808 1887 1663 1738 1835 1928 1667 1768 1891 
1666 1739 1813 1888 1671 1743 1840 1929 1668 1770 1892 
1670 1742 1815 1889 1672 1744 1841 1930 1671 1773 1907 
1673 1745 1819 1906 1678 1748 1845 1932 1672 1784 1913 
1674 1746 1820 1909 1679 1752 1847 1933 1678 1787 1920 
1675 1750 1821 1912 1681 1768 1854 1935 1679 1798 1924 
1679 1753 1826 1914 1683 1771 | 1855 1939 1681 1799 1926 
1680 1756 1831 1916 1686 1774 1862 1942 1686 1814 1928 
1684 1757 1832 1921 1687 1775 1866 1944 1687 1823 1929- 
1685 1761 1839 1926 1690 1778 1877 1948 1690 1834 1933. 
1696 1763 1841 1931 1693 1784 1879 1952 1691 1835 1935 
1698 1767 1845 1935 1696 1785 1889 1953 1693 1836 1939 
1701 1771 1847 1938 1703 1786 1891 1956 1699 1837 1942 
1705 1773 1849 1943 1707 1796 1893 1959 1706 1847 1944 
1709 1778 1851 1945 1710 1802 1894 1962 1709 1854 1946 


1711 | 1781 | 1855 | 1949 | 1711 | 1804 | 1901 1711 | 1856 | 1952 
1715 | 1782 | 1859 | 1950 | 1712 | 1805 | 1907 1716 | 1857 | 1953 
1720 | 1786 | 1860 | 1951 | 1713 | 1807 | 1913 1723 | 1858 | 1955 
1721 | 1787 | 1866 | 1954 | 1714 | 1809 | 1917 1725 | 1861 | 1959 
1725 | 1794 | 1868 | 1956 | 1716 | 1810 | 1918 1735 | 1862 
1728 | 1797 | 1870 | 1958 | 1718 | 1814 | 1920 1739 | 1868 


1729 1802 1873 1959 1722 1818 1922 1744 1876 
1731 1803 1878 1725 1832 1924 1752 1877 
1733 1806 1883 1737 1833 1927 1755 1885 





据 此 ,以 1663—1962 年 资料 分 别 统计 水 , 旱 , 蝗 的 发 生 次 数 ,用 列 联 表 进 行 分 析 ( 表 2、 
3、4、5)。 FEES LURIA, TERA 1) 水 时 年 与 一 般 年 之 间 的 飞 蝗 发 生 数 次 有 显著 差异 (2 
式 ); 2) 时 年 对 飞 蝗 大 发 生 的 作用 又 远 较 洲 年 或 讲 后 第 二 年 为 大 (3 式 、4 式 ); 3) 游 后 


R2 不 同 降水 年 份 与 蝗虫 发 生 关系 比较 RI 水 .早年 与 昕 虫 发 生 关系 比较 
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第 二 年 与 涝 年 比较 , 则 前 者 的 作用 又 较 后 者 显著 (5 式 )。 


Rng \2 
3 2 ( Rik 一 at) 
= 300 >) » 一 一 


i=] k=1° nink 
P(X? > 29.91) < 0.001 


= 29.91 (2) 


X= n( ayn 一 nanny = 23.9 (3) 
NNN. n.2 
PCH? > 23.9) < 0.001 X3 w = 10.82 


12 = Pinna — nvr)” g4 (4) 

i 力 1.722.77.112 .2 

y = Omna — nmm) _ 4 6 (5) 
y.22.N.4N.2 


P < 0.05 P(0.05) = 3.8. 

又 用 近 50 年 (1913 一 1962) 浴 河流 域 各 月 降水 等 级 图 的 资料 ? 将 各 季 的 日 讲 情况 ?3 与 
同期 蝗虫 发 生 资料 结合 进行 概率 分 析 ( 表 6)， 进 一 步 说 明 夏 季 (6 一 8 AFB CAR 
生 的 影响 最 大 ,吻合 条 件 构 率 达 0.846, 春季 次 之 ,秋季 与 冬季 的 干旱 对 飞 蝗 影 响 不 大 ， 而 
春季 中 又 以 5 月 干旱 与 飞 蝗 发 生 的 条 件 概率 最 高 。 


R6 不同 季节 旱情 与 飞 蝗 发 生 次 数 比 较 






发 生 干旱 的 季节 


春季 (3 一 5 月 ) 旱 年 为 10 次 

夏季 (6 一 8 月 ) 旱 年 为 13 次 
秋季 (9 一 11 月 ) 旱 年 为 7 次 
冬季 (12 一 2 月 ) 旱 年 为 11 次 







(te 50 年 中 无 条 件 概率 了 = 6 一 0.32 ) 


早 与 涝 都 是 降水 及 温度 等 复合 因素 的 指标 , 须 分 别 找 出 其 作用 的 范围 或 数值 ,以 5 一 
6 月 及 7 一 8 月 降水 量 分 别 与 夏秋 蝗 发 生 级 数 进行 分 析 , 明确 夏 蝗 大 发 生 的 降水 条 件 系 
5 一 6 月 降水 量 在 100 毫米 以 下 , 秋 蝗 的 指标 系 7 一 8 月 降水 量 在 210 EKAT, 以 列 联 表 
进行 显著 性 检验 ,二 省 结 果 均 极 显著 , 分 别 为 如 一 11，P( 妇 > 11) < 0.001, X = 7.693 
P(X > 7.693) < 0.02。 再 以 各 年 5 一 6 月 及 7 一 8 月 的 平均 温度 代表 该 时 期 的 温度 条 件 ， 
UFR (z = Ke) 代表 该 时 期 的 湿度 条 件 ， 结 合 蝗 情 进行 图 析 (结果 见 图 1) 得 
知 ,夏秋 蝗 大 发 生 的 温 湿 庆 条件 是 干 湿 系 数 < 3 并 且 < y 一 0.95cotzr*， 而 一 般 发 生年 的 
于 湿 系 数 是 > 3M > y = 0.95e%", . , 


1) 和 月 降水 等 级 图 资料 条 由 中 央 气 条 研 究 所 供给 ,该 图 将 降水 分 为 i 1 级 为 大 水 ,2 PARE 3 级 为 中 常 ， 
4 RAR, 5 级 为 大 旱 。 

2) 我 们 确定 该 区 每 季 中 有 一 个 月 的 降水 等 级 指标 为 5 ,而 其 它 .2 月 均 在 3 级 以 上 者 为 旱季 ,每 季 中 有 一 个 月 为 1 
级 而 其 它 2 个 月 在 3 以 下 者 为 水 季 。 
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Cy 二 .蝗虫 发 生 期 间 的 月 平均 温度 ，x 一 蝗虫 发 生 期 间 的 月 干 温 系 数 。) 





R/S 
图 1 suum TR 
+ 大 发 生年 一 发 生 轻 度 年 

夏 、 秋 凰 发生 级 数 及 发 生 数 景 除 与 上 述 旱 渤 程度 有 关外 ， 倚 有 其 它 环 境 因 素 的 作用 ， 
Al AMEE AREA, A 1950—1963 年 5 一 9 月 份 的 月 平均 温度 、 相 对 湿度 、 降 水 量 、 
降水 日 数 、 湖 水 位 、 上 代 虫 口 、1 月 的 温度 、 越 冬 卵 及 其 死亡 紊 分 别 对 东亚 飞 蝗 的 年 发 生 
级 数 ,夏秋 凰 发 生 级 数 与 发 生 数 晤 以 及 各 因素 之 间 进 行 了 相关 距 阵 的 分 析 ， 所 得 结果 如 
下 : 

(1) 夏 蝗 发 生 级 数 与 1 月 平均 温度 呈正 相关 〈?* 一 0.55) ， 与 闭 水 位 呈 负 相关 C= 
一 0.64), 与 越冬 卵 量 级 数 呈 正 相 关 (r = 0.64),5 6 月 降水 呈 负 相关 (7 = 一 0.55), 以 上 > 
值 均 在 P <0.05 水 平 。 秋 蝗 发 生 级 数 与 潮水 位 成 负 相 关 〈* = 一 0.53)， 与 7 HKE 
负 相 关 (> = 一 0.52), 显著 性 测定 P < 0.05; 越 冬 卵 死亡 众 与 9 月 降水 量 呈 正 相 关 (7r 一 
0.802), 而 越冬 卵 量 与 湖水 位 晨 负 相 关 (r = 一 0.746), 显 落 性 测定 P< 0.01。 

(2) BARAZER SRAI RIE (r = 0.811), 51 月 温度 本 呈正 相关 (+ 一 
0.55)， 潮 水 位 与 越冬 卵 景 呈 负 相关 (r = 一 0.785)，5、6 月 降水 日 数 与 湖水 位 星 正 相关 
Cr 一 0.75)。( 显 落 性 测定 除 1 月 温度 在 P < 0.05 外 ,其 它 均 在 P < 0.01.) 

(3) 夏 蝗 在 湖 认 上 的 发 生 面 积 与 4 月 或 5 月 的 湖水 位 呈 极 显著 的 负 相 关 (r 二 一 0.72) 
Cr = 一 0.79), 秋 蝗 在 湖 滩 上 的 发 生 面 积 , 则 与 7 月 的 平均 湖水 位 星 极 显著 的 负 相 关 〈* 一 
一 0.87)o 

整个 湖区 的 夏 蝗 发 生 面 积 与 前 年 9 月 及 当年 4 月 平均 湖水 位 及 当年 4 月 与 5 月 最 高 
湖水 位 ， 或 前 年 9 月 最 高 水 位 及 10 月 平均 湖水 位 的 相关 关系 均 在 P > 0.05， 相 关 不 显 
落 。 而 秋 星 发 生 面 积 则 与 7 月 至 8 月 降水 量 呈 负 相 关 , 显 著 性 检验 P < 0.05, 

(4) 本 区 各 环境 因素 之 间 的 相关 则 有 : 5 月 温度 与 8 月 降水 呈 负 相关 ,5 月 降水 与 6 
月 降水 呈正 相关 ,5 月 温度 与 8 月 温度 呈正 相关 ，5 月 温度 与 6 月 温度 呈正 相关 ，6 月 温 
度 与 7 月 温度 呈正 相关 , 7 月 温度 与 7 月 降水 呈 负 相关 。 

种 从 生物 学 因素 的 分 析 

种 恒生 物 学 因素 包括 两 方面 , 朗 种 全 动态 的 时 间 ( 或 历史 的 ) 序 烈 及 种 全 密度 因素 ,前 
者 表现 于 前 后 若干 代 或 更 长 时 间 的 持续 影响 ， 即 种 下 数量 的 增 减 是 否 有 前 后 延续 的 周期 
性 波 动 关系 ,后 者 表现 于 因 密 度 增 减 所 产生 的 种 轻生 殖 力 及 成 活 率 的 差异 , 序 当 代 密 度 与 
TREER AMT: 

1) 东亚 飞 蝗 的 发 生 周 期 性 问题 
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RW ES ON FAES RNA ME, KRAP E ey 
te PEK Aes —, WEES ORAM RT ES MS hese 
REMEE BA OT, AAA ABE (DER 1956、1958) , 找 出 大 发 生 时 间距 离 最 大 
间隔 为 9 一 11 年 ,最 小 间隔 为 4 一 5 年 ,在 二 次 大 发 生 之 闻 如 时 间距 痪 多 于 5 年 , 则 其 则 有 
1 一 3 次 小 发 生 ， 每 次 小 发 生 距 议 为 2 一 3 年 。 为 了 进一步 证 明 此 等 时 间 间 隔 是 否 有 规律 
性 ,以 便 在 长 期 预测 中 运用 ,又 将 1913 一 1962 年 的 50 年 飞 蝗 发 生 资料 进行 了 谐 波 分 析 及 
积分 曲线 分 析 ， 根 据 谐 波 分 析 结 果 ( 图 2) 可 以 看 出 东亚 飞 蝗 的 发 生 周 期 有 2 一 3 年 , 5 一 6 
年 及 25 年 左右 的 三 个 周期 。 


1.302 
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0.807 
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1 2 5 6 7 8 9 10 11 12 1314 15 16 17 19 20 22 23 
波 F 
图 2 东亚 飞 蝗 种 学 动态 的 谐 波 分 析 
”再 则 根据 积分 曲线 分 析 结果 , 似 亦 存在 有 25 年 的 长 周期 (图 3)， 而 此 周期 又 与 本 区 
5 一 8 月 旱 月 的 周期 相 购 合 (图 3)。 
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1913 1917 1921 1925 1929” 1933 1937 1941 195 199 1953 1957 . 1961 


年 a 
图 3 ， 东 亚 飞 蚊 利 生动 态 的 积分 曲线 分 析 及 其 与 水 蝗 区 5 一 8 RE : 
HEME 5 一 8 月 旱 月 积分 沿线 。 一 一 东亚 飞 妈 的 积分 曲线 


利用 积分 曲线 讨论 周期 间 题 ,在 气象 方面 已 有 先例 (Tapc, 1956; 杨 鉴 初 ， 1958)， 以 
此 用 以 分 析 长 周期 的 结果 亦 谷 满意 ， 措 且 积 分 曲线 与 谐 波 分 析 在 长 周期 方面 所 表现 的 结 
果 亦 大 致 吻合 ,只 是 谐 波 分 析 的 波 序 在 长 周期 方面 年 距 较 长 而 且 在 20 年 以 上 的 周期 仅 有 
25 年 一 个 。 因 此 在 进行 此 项 分 析 时 ,采用 两 种 方法 以 互补 不 足 。 

从 上 还 结果 可 以 看 出 东亚 飞 幅 由 于 具有 二 个 以 上 的 周期， 各 周期 间 互 相交 错 , 因此 ， 
多 观 分 析 在 发 生 级 别 上 就 表现 出 不 明显 的 周期 。 

2) 种 堆 的 密度 因素 
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东亚 飞 蝗 不 同 密度 下 的 生殖 力 : 以 夏 , 秋 蝗 不 同 密 度 种 从 生殖 力 的 结果 与 Watt(1960) 
的 高 维 超越 曲线 方程 进行 配合 ,完全 吻合 , 夏 蝗 生殖 力 高 于 秋 凰 ， 但 秋 蝗 出 现 生殖 力 高 党 
的 理论 值 的 密度 则 低 于 夏 蝗 ( 另 广 讨论 )o 

种 者 死亡 率 与 密度 的 关系 ,根据 卵 与 蝗 晴 分析 的 结果 , 则 分 别 吻合 于 高 次 函数 曲线 方 
程 与 超越 曲线 方程 ( 另 交 讨论 )。 ' | 

此 和 外 ;从 上 述 大 发 生 之 间 的 独立 性 测定 结果 中 ,区 已 说 明 种 价 数 量变 化 的 前 后 持续 联 
Ro | 

人 为 因素 的 作用 分 析 

3) 人 为 因素 的 作用 包括 两 方面 , 郎 化 学 防治 和 土地 利用 ， 前 者 系 影 响 虫 口 基 数 的 决 
定性 因素 ,后 者 对 发 生 密 度 及 发 生 面积 都 超重 要 作用 ， 由 于 化 学 防治 的 结果 ,可 以 虫口 基 
Beal BS hs (A oe) Se BERS, 士 地 利用 程度 与 面积 则 直接 决 定 于 前 水 位 及 内 游 面 积 ,因此 
交接 采用 湖水 位 及 内 涝 面积 的 数值 较 用 土地 利用 级 数 更 为 方便 。 


二 、 几 种 预测 法 的 建立 


从 上 述 对 东亚 飞 蝗 作用 因素 的 分 析 结 果 ; 结 合 东 亚 飞 蝗 种 移动 态 的 特征 ,作者 在 建立 
东亚 飞 蝗 的 预测 式 时 ,从 下 列 三 方面 进行 考虑 : 

1) 根据 东亚 飞 蝗 在 历史 上 种 司 变 动 的 时 间 序 列 的 特征 ， 建 立 预 测 模 式 ,这 是 考虑 到 
种 从 在 长 期 内 一 么 列 连 续 的 发 生 数 量 或 强度 ,全 面 地 反映 出 该 种 谷 在 历史 上 的 变化 ,这 种 
变化 包 合 一 切 内 在 与 外 在 因素 的 作用 总 和 ,因此 可 以 利用 此 规律 进行 预测 。 

2) 根据 种 众生 物 学 特性 的 长 短 周 期 规律 结合 随机 影响 的 作用 ,来 建立 预测 方程 。 

3) 根据 大 量 外 界 作 用 因 来 对 种 谷 变 动 的 影响 建立 预测 方程 。 

因此 ,上 述 三 类 预测 式 代 表 三 个 不 同 的 方面 ,并且 三 种 方法 各 有 其 长 , 亦 各 有 不 足 , 员 
为 三 个 独立 之 预测 模式 ， 但 应 将 其 视 为 组 成 总 预测 式 的 三 个 相互 依赖 的 部 分 。 兹 分 述 如 
TF: 

(一 ) 根据 历史 规律 性 的 分 析 进 行 预 测 

以 前 述 种 恒 各 年 间 的 发 生 级 数 为 依据 , 绘 成 发 生 级 数 变 化 序列 图 (图 4), 图 中 升降 序 
列 趋势 代 表 种 从 在 综合 因素 作用 下 的 动态 型 (Population dynamic form), 此 法 系 根据 动态 
型 的 规律 性 趋势 ,进行 外 推 估 值 。 历 史 规 律 性 特征 系 从 下 烈 几 方面 进行 分 析 与 处 理 。 





0 
“1913 1917 1921 1925 1929 1933 1937 1941 “1945 1949 1953 -1957 1961 


年 份 
图 4 1913—1961 RY CER A Ab FA hes HR 


1. RARER: 为 了 集中 起 见 把 1 RHA 2 级 , 10 RHA 9 级 ， 然 后 
利用 马 氏 链 的 一 步 转移 矩阵 ,把 50 年 的 资料 进行 分 析 , 结 果 如 Pij 
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0.14 0.50 0.07 0.14 0.14 0 0 0 
0.40 0.10 0.10 010 10 0.20 0 0 

Py Py +++ Py 0.1250.25 0 °0 0.1250.25 0.25 0 

p, Pa Pas a Ps» _| 033 0 033 0 0 0.33 0 0 
: : -0.60 0 040 0 0 0 0 0 

Po Pa- Po 0.33 0.17 0.33 0 017 0 0 0 

0 0 050 0 0 0 0 050 

100 0 0 0 0 0 0 0 


FREA TERETI AP ARREN, B 表示 下 降 , C 表示 平稳 , 则 其 转移 
SE RHE: 
0.86 0 0.14 
0.50 0.40 0.10 


P24 Pag Pac 0.55 0.45 0 

P34 Pag Pac 0.33 0.67 0 

Pu=) >: 217] 020 oso o 
Pos Pox Poc 0 1.00 0 

, 0.50 0.50 0 

0 100 0 


从 第 一 个 矩阵 的 结果 ,可 以 看 出 2 级 转 入 3 级 的 概率 最 大 , 4 级 转 人 3 级 、7 级 、8 级 
的 概 牵 较 多 ,第 二 矩阵 的 结果 则 显示 出 2 级 上 升 的 概率 最 多 , 5 级 以 上 至 7 SY PARE 
BRD o 

2. 发 生 状态 的 持续 波动 分 析 将 50 年 的 蝗虫 发 生 级 数 进行 平均 , 求 进 平均 级 数 (7)== 
4.14, 然 后 从 图 4 中 可 以 看 出 下 烈 趋 势 。 

1) 当 级 数 超 过 距 平 土 2 时 , 则 原来 的 升降 区 势必 然 发 生 逆转 , HY SBI 96 % 0 

2) Ea r 表示 第 i 年 的 发 生 级 数 , 从 图 4 可 知 ;, 若 ria tia y 则 wa < tis 
#2, > xia > xia, W r > x20 #EE -REAT ras < 4.1， 以 上 吻合 概 牵 为 
82% . 

3) 若 xS 3, Xitl < 3, 则 Xian 之 Xip 吻合 概率 为 89%。 

3. 发 生 型 特征 的 相似 性 分 析 

在 长 期 的 波动 趋势 中 ,存在 有 某 几 年 的 波动 内 型 与 另 些 年 相似 的 情况 ,如 1916 一 1930 
年 之 间 及 1949 一 1962 年 之 间 的 崇 型 等 是 。1916 一 1921 年 、1922 一 1927 年 及 1928—1930 
4p HUA A IESE 2 次 后 序 行 下 降 , 而 后 出 现 一 单 罕 ; 1949—1956 年 与 1957—1962 
年 的 党 型 系 前 峙 连续 上 升 3 次 及 下 降 2 次 后 出 现 一 单 罕 。 

KEMAH 1963 年 的 发 生 级 数 : 根据 转移 和 矩阵 的 分 析 ， 共 已 知 1962 年 发 
生 级 数 系 3 级 ,从 前 一 转移 和 矩阵 中 ,可 以 看 出 3 级 转 入 2 级 的 概率 最 大 。 根 据 罕 型 相似 性 
分 析 , 1957—1962 4 HAIG 1950—1956 年 罕 型 在 发 展 上 基本 一 致 ,好 前 案 后 应 出 现 一 单 
#1 1963 年 正 是 单 罕 的 终点 , 因此 , 估计 1963 年 的 发 生 级 数 应 系 1 一 2 级 ,属于 小 发 生 
年 。 . . 
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(=) 应 用 随机 序列 与 周期 方程 的 预测 

以 1913 一 1962 年 的 50 年 蝗虫 发 生 级 数 资 料 , 按 年 序 组 成 一 时 间 序 列 (y,，: 二 L, 2, 
… ,50) BPR CAA 25 年 以 上 及 10 年 以 下 的 长 短 周期 ， 因 此 假设 此 序列 随时 间 
变化 的 因素 ,和 柔 由 三 类 影响 所 组 成 ; 1) 长 周期 影响 , 2) 短 周期 影响 (10 年 以 下 的 周期 )， 
3) 随机 影响 ,如 此 ,其 表达 式 朗 可 写成 

yi =F) + ECE) + gle) (6) 

AR yii th ee Be OA EEL, gO SAH, 

从 累积 的 蝗虫 发 生 资 料 中 将 长 周期 项 (趋势 项 ) 及 随机 项 逐步 分 离 ， 再 结合 周期 项 的 
变化 规律 进行 分 析 后 作出 外 推 。 

1. 长 周期 项 的 分 高” 采用 移动 平均 法 消除 与 减弱 较 短 的 周期 波动 ,按照 24 十 1 个 的 
相 邻 项 求 其 平均 





/ 1 
Yau = R41 Cy. ot ++ + vanes) (7) 


其 中 yaar 是 2k 十 1 年 的 移动 平均 值 ,其 坐标 沙 于 十 1 年， 以 不 同 & 值 试 作 的 结果 ,发 现 
X k = 6 k, Èp 13 年 移动 平均 ,趋势 比较 明显 (图 5 )。 此 曲线 可 以 用 

JA) = a+ bx csin(wi +a) (8) 
方程 来 逼近 。 





20 25 30 35 40 45 50 
年 数 
图 5 ARIE ese 
一 一 和 一 6 移动 平均 曲线 — FO) 曲线 


即使 满足 J(a, b,c) 一 5 {z,—[a+ bx + csin (wt + a)]} 为 最 小 , Hh z, B y 
1=7 


的 13 年 移动 平均 值 。 


Jla, b, c) = 0 
ða | 


OJ(a, b, c) = 0 
Ob 


Jla, b, c) = 9 
Oc 


解 得 a = 4.475, b = — 0.0084, c = 0.38, Hh w, a 由 图 中 直接 可 得 w= Za =0 
代入 (8) 式 即 得 


> 
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f (2) = 4.475 一 0.00842 + 0.38 sin (Z ‘) (9) 


”2. 随机 项 分 离 ST FEL Be FAMERS (Araon, 1952) 来 处 理 随机 影响 。 随机 项 
CE) = y — C) o l (10) 
BD AAAS BRAT HD BS HH i BJ FF 
把 &, 作为 一 平稳 随机 序 烈 ， 进行 线性 外 推 处 理 ， BN ¢ 年 的 随机 影响 与 其 前 《年 的 
随机 影响 成 线性 画 数 关系 。 


E = aber + obi Fee + Ane pk (11) 
IRERE R, REN = 2 BEE, EN a 
EF = ad- 十 oo- (12) 


Bam) = E| 6 一 Sue —s) Geke gaga) 


iTi 
= Bo 20Bey — 2a:B m + Bo 十 2aoxBa + Bw) 
为 使 fC ayo) 还 到 极 小 ， 满足 条 件 





Of (ara) =o. 
Oa 
Of Caa) _ 
Oa E 
得 
aB + aBa = -| (13) 
aBa + Bo = Bo 
在 方程 (13) 中 
Bo 一 二 >ë 
n tol 
1 5 
Bo = —— >) ëm (14) 
1 一 工 ?21 
Bw 一 EEn J 
n— 2z 
性 处 > 一 50, 由 方程 (13) 解 得 
a = —0.1186, æ, = —0.1452 
TAG AMX 
E*(2) = —0.1186F(2 — 1) — 0.1452E(2 — 2) (15) 
3. 周期 项 的 处 理 M y 中 分 离 出 趋势 项 及 随机 项 之 后 ,其 余 项 用 短 周期 项 处 理 即 
ye — fC) 一 全 (一 8 (16) 


由 于 在 趋势 项 中 f(z) 作 了 13 年 滑动 平均 , NUL APM EeSMPoAs 
出 去 ， 现 在 考虑 的 应 为 2 一 12 年 , 并 由 于 2 一 6 年 分 别 包含 在 8、9、10、12 的 周期 中 综合 
考虑 ,所 以 仅 需 对 7 一 12 年 周期 影响 加 以 处 理 , 邹 对 以 下 序列 进行 分 宛 。 

g° (2) = gr@) + gG) + +++ + goe) (17) 
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(18) 


gay = fP e a) bf) xt) A xt + xf) 
gh = xy) 十 x 十 x 办 十 x 从 十 x 人 十 ar 


ee 


Jt 50 个 方程 (最 后 2 年 留 作 外 推 检验 ,实际 用 48 个 ,其 中 
zx 多.… x 为 7 年 周期 全 部 画 数 值 
te «++ x 为 8 年 周期 全 部 图 数 值 
依次 类 推 ,由 于 (18) 式 中 , 便 有 11 个 未 知 数 是 不 独立 的 ,如 
7 年 周期 中 有 1 个 不 独立 画 数 
8 年 周期 中 有 +4 个 不 独立 画 数 
9 年 周期 中 有 3 个 不 独立 画 数 
10 年 周期 中 有 2 个 不 独立 男 数 
11 年 周期 中 有 1 个 不 独立 画 数 


xh =0 


因此 合 


1 = P = @ = ) = 9 


xD = x) = x2) = 0 
1) = xP = 0 
xD = 0 


TES — AE BaP = P = h l xf) = PP ABM, TEE 47 个 方程 中 只 有 45 个 
未 知 数 , 郎 成 为 解 了 矛盾 方程 问题 ,这 样 将 上 上述 方 程 变 为 


g2 = xy) cb ext? 十 x? 
ga = xP + x8 + 2 + x 
g = xs? + xf) + xf) + af) + xf? 
gs =) + xe) + eo) + xe + xf + xO (19) 


ge = af + OF A + AO + a + 2 





1917 1922 1927 1932 1937 1942 1947 1952 1957 1962 
. 年 份 
图 6 随机 序列 及 周期 的 预测 值 对 东亚 飞 蜂 实际 发 生 值 的 内 插 外 推 比较 
一 :一 实际 发 生 值 一 一 六 十 人 十 8 的 理论 估 值 ---- 外 推 值 
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《19) 式 在 电子 计算 机 上 用 迭代 法 去 解 , 朗 可 求 得 此 等 未 知 周 期 的 近似 值 。 
上 述 各 项 的 年 值 及 其 和 (部 预测 值 ) 详 见 表 7 所 列 结果 。 
从 表 7 结果 看 出 , (I) 根据 y 一 f OHE O He” O 方程 对 48 年 资料 进行 内 播 的 结 
果 , 可 以 清晰 地 大 出 理论 值 与 实际 佳之 间 的 吻合 程度 良好 。(I1) 根据 y = FG) 十 EO 
R7 不 同 处 理 与 实况 比较 





A | | 随机 项 | AAR |E | 1G) +E) +22) 


实 

1 ` ; 

年 别 | (安生 级 数 ) | fC?) EC) a(t) 
1913 7 4.63 0.269 2.323 4.9 7.3 
1914 4 4.53 0.649 | 一 0.883 5.18 4.3 
1915 3 4.42 0.716 | -1715 5.14 3.4 
1916 1 4.31 0.439 | 一 3.166 4.74 1.6 
1917 2 4.23 0.447 | 一 1.865 4.67 2.8 
1918 3 4.2 —0.028 | —0.413 4.17 3.8 
1919 2 4.1 0.784 | —0.739 | 4.78 4.0 
1920 6 4.09 0.423 2.364 4.51 6.9 
1921 2 4.1 0.078 | 一 1.244 4.18 2.9 
1922 3 4.12 —0.028 | —0.269 4°09 3.8 
1923 4 4.18 0.438 0,156 4.62 4.8 
1924 6 4.25 0.184 1.741 4.43 6.2 
1925 2 4,33 一 0.182 | —1.987 4.15 2.2 
1926 5 4.42 0.022 0.117 4.44 4.6 
1927 4 4.52 0.270 | —1.398 4.79 3.3 
1928 8 4.60 —0.023 2.422 4,58 7 
1929 0 -4.67 —0.328 4.592 4,34 9 
1930 2 4.73 —1.226 | 一 2.947 3.604 0.7 
1931 2 4.75 一 0.450 | —3.535 4.13 0.8 
1932 4 4.75 0.723 | 一 2.822 5.47 2.7 
1933 7 4.73 0.488 2.716 5.21 8 
1934 4 4.69 —0.398 | —1.124 4.29 3.2 
1935 8 4.62 一 0.538 3.353 4.08 7.5 
1936 4 4.53 -0.301 | -0.573 4.23 3.7 
1937 2 4.43 —0.428 | 一 1.975 4.00 2 
1938 3 4.32 0.365 | 一 1.433 4.69 3.2 
1939 7 4.22 0.509 2.628 4.73 7.4 
1940 2 4.11 —0.138 1.483 3.97 5.4 
1941 3 4.02 —0.153 | —0.466 3.87 3.4 
1942 6 3.95 0.427 2.127 4.38 6.5 
1943 2 3.91 —0.095 | —1.730 3.82 2.1 
1944 5 3.89 —0.071 1.484 3.82 5.3 
1945 2 3.90 0.146 | —2.096 4.05 2 
1946 6 3.92 0.064 1.878 3.98 5.9 
1947 4 3.97 0.029 | —0.263 4.0 3.7 
1948 3 4.05 —0.306 | —0.848 3.74 2.9 
1949 2 4.13 0.120 | 一 2.692 4.25 1.6 
1950 3 4.22 0.287 | —1.686 4.51 2.8 
1951 5 4.31 0.309 0.266 4.62 4.9 
1952 7 4.40 0.095 “| 2.435 4.5 6.9 
1953 6 4.47 一 0.409 1.795 4.06 5.9 
1954 4 4.53 一 0.571 0.074 3.96 4.1 
1955 7 4.55 一 0.174 2.634 4.38 7 
1956 2 4,55 —0.214 | -2.241 4.34 2.1 
1957 3 4.53 —0.053 | —1.362 4.48 3.1 
1958 3 4,49 0.552 —1,892 5.04 3.1 
1959 7 4.42 0.399 2.272 4.82 7.1 
1960 3 4.33 —0.09 — 1.025 4.24 3.2 
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+ g#( 人 方程 ,用 1913 一 1960 年 资料 、1913 一 1961 年 资料 分 别 对 1961 年 、1962 年 发 生 情 
驶 进行 外 推 检 验 , 结 果 亦 比较 满意 (图 6 )。 

(=) SRRBF ARMS 

在 上 述 ( 一 )( 二 ) 两 种 方法 中 , 均 柔 由 历史 资料 直接 进行 预测 ,不 考虑 任何 具体 变量 的 
作用 ,实际 上 ,自然 界 中 有 许多 对 蝗虫 数量 起 作用 的 环境 因素 ， 它 们 的 变化 及 作用 具有 时 
间 及 空间 特征 ,因此 在 应 用 上 述 两 汶 法 的 同时 ,还 应 考虑 此 类 环境 因素 的 具体 变化 对 蝗虫 
数量 (或 级 数 ) 所 起 的 作用 ,作者 在 计算 所 电子 计算 机 的 帮助 下 , 应 用 多 因素 逐步 直线 通 归 
与 曲线 迎 归 分 别 进行 了 东亚 飞 蝗 数值 (或 级 数 ) 的 预测 与 内 插 的 比较 ， 从 而 确定 出 一 种 在 
这 方面 较 好 的 预测 方法 , 兹 分 述 如 下 . 

1. 多 因素 深 步 直线 通 妇 预测 法 (1963—1964 年 外 推 各 值 系 借用 史 久 因 、 惟 玉 簿 编 的 
程序 。) . 

在 给 出 的 大 量变 量 中 (此 次 给 出 的 最 多 变量 是 14 个 ), 通 过 给 定 的 信 度 检验 ， 权 衡 关 
示 大 小 ， 将 其 中 与 星 虫 数量 (或 级 数 ) 关 系 较 小 或 可 用 另 一 些 因素 值 代 蔡 的 茶 些 因 来 值 滤 
法 ,以 留 下 的 因素 值 构成 凶 归 式 方程 


y = by + Bri bixa + °° + Ox, (20) 
7 一 6 十 pixi CF RH) 


?9 一 bo + bixi + byx2 (F 检验 ) (21) 


EE FA LVR PRA, A SER, PE 1950—1963 5 1952—1963 年 资 
料 , 分 别 用 级 数 与 密度 进行 处 理 , 双 每 种 处 理 均 为 二 项 (1) 包 括 1963 年 资料 外 推 1964 年 ， 
《2) 不 包括 1963 年 资料 外 推 1963 年 。 

















所 得 各 项 预测 式 如 下 : 
1) SERRA Be 
(1) 外 推 1963 年 EKF F = 1.5 (F = 2.5 下 的 公式 相同 ) 
y = 43.24 一 2.03xs 十 4.953xy (22a) 
(2) 外 推 1964 年 fake F = 1.5 TF 
y = 287.8 一 0.273x, 一 33.9x5 + 3.703x7 + 9.254xo (22b) 
HEKE F = 2.5 T 
y = 14.05 + 4.91 x; (22c) 
式 中 y—— AWE, aó HEKE, aDC, r6 月 湖水 位 ，x: 一 一 越 
AHE, xo 5.6 月 降水 日 
2) 夏 蝗 发 生 级 数 


(1) 外 推 1963 年 下 一 1.5 
y = 10.68 十 0.144xi 一 0.812x? 十 0.309xz5 一 0.002x6 一 0.016xy 十 0.276xs 一 0.383xe (23a) 
F=2.5 下 
y = 10.797 + 0.114%, 一 0.973xa 十 0.274xs (23b) 


1 Ame, x 湖水 位 ，xs 一 一 5 月 降水 日 ， Xe 6 月 降水 日 ， X7 5 月 














RFR xy 
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降水 量 ，xs 一 一 6 月 平均 温度 ，xs 一 5 月 平均 温度 。 
(2) 外 推 1964 年 下 一 15 (F 二 2.5 同 ) 
y 一 一 4.59 一 0.945x, 十 0.456x3 + 0.482x4 十 0.386x5.— 0.009xy 十 0.522x5 . (23c) 
REP x WAZKAL, x3 上 代 虫 口 级 数 ，z 一 一 6 月 降水 日 数 ，xs- 一 5 月 降水 日 数 ， 
5 月 降水 量 ，xs 一 一 6 月 降水 量 。 
3) KERERE 
(1) 外 推 1963 年 F=15 (F = 2.5 Ñ) 

















7 




















y = 806.8 一 0.752x, 一 54.582x3. (24a) 
TAA x4 湖水 位 ，xs 一 一 8 月 降水 日 数 。 i 
(2) 外 推 1964 年 F= 1.5 
y = 一 8.523 一 0.451x, 十 0.309x, 一 49.25xa 十 23.278x6 (24b》 
SAR. xy 8 BEZK, x 7 月 降水 ， x3 HAKEE,- xs 一 一 7 月 平均 温度 。 


4) 秋 蝗 发 生 级 别 
(1) 外 推 1963 年 下 一 1.5 








y = 一 9.062 一 0.415x, + 0.632xy 一 0.006x6 (25a) 
F =2.5 
y. = 一 15.523 — 0.007 x; + 0.69x, (25b>} 
AP xy 潮水 位 ，x 7 月 降水 ，xs 一 一 8 月 平均 温度 。 


(2) 外 推 1964 年 下 一 1.5 
y 一 32.55 一 1.555x1 一 1.78x, + 0.159x3 一 0.375x4 — 0.008x6 



































十 1.05x, 一 1.088xs (25c) 
F=25 
y = 一 14.7 — 0.006x, + 0.6512, (25d) 
式 中 xy WAZKAL, x2 ERRORS, x3 8 AMKA, x, 7 月 降水 日 , xs 
8 BKE, re 7 BKE, x 8 月 温度 ，xs 7 月 温度 。 
5) 年 发 生 级 数 
(1) 外 推 1963 年 
y = 一 26.89 一 0.742x; 十 0.326x6 + 0.564xu + 0.85973 (26a) 
(2) 外 推 1964 年 
y 一 一 30.01 一 0.015x,; 十 0.639xu 十 0.708%x43 (26b) 
式 中 ao BEKE, aiK, xs 一 一 5 月 降水 日 ,za 一 一 8 月 平均 温度 ,xu 一 一 
7 月 平均 温度 。 | 


NEESER 1963 年 的 外 推 与 内 插 结果 见 表 8 及 表 9 。 
2. FARA th Bie A 
公式 
y = ag t aixl 十 a2x2 + +++ + a,x, 
二 bb (27) 


-十 cixixa 十 cHxlxa 十 十 CjjxXixj 
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R8 多 元 直线 方程 内 插 喘 莽 比 较 
级 数 RK Hl 测 
年 别 年 发 生 级 数 Kk em 级 数 tk g 


F=1.5 | F=2.5 | F=1.5 


1950 .075 1.238 一 0.421 . 。 一 


1951 .061 0.426 | —0.071 . . o . | 一 

1952 .075 2.794 | 一 0.197 . . . . 9.77 
1953 .012 | -0. 1.126 . . . . 104.71 
1954 .008 0.026 .048 .122 | . . . 0 

1955 .046 1.295 .717 . . —14, .82 一 3.24 
1956 | —0.003 .296 .089 . . . . 2.79 
1957 .145 . . .1 . . . 8.00 
1958 . . .934 . .25 . .835 22.49 
1959 . . .357 . . . . 36.35 


1960 。 ` -499 - . - .505 一 43.89 
1961 . 。 .754 。 . 。 . — 126.25 
1962 . 0. . -828 。 ` 。 。 一 28.34 








级 «BY i 数 值 
年 发 人生 级 数 Ze 级 数 KERK E 星 fk 


1 通 归 预测 值 1 RIE E WIS Mw ie 1783 WAH 1763 通 归 预测 值 
SEL | F=1.5|F 一 2.5| 实 驶 |F=1.5|F 一 2.5| 实 现 | F=1.5| P= 2.5 | RI | F=1.5 | P= 2.5] K| F=1.5 | P= 2.5 






y= aot aay t+ biri CF 检验 ) 
y= ag + ax + aze + biri + bt + cear, 《检验 ) (28) 


夏 蝗 数值 外 推 63 年 .下 一 1.5 


y 一 一 0.125 一 0.009%, + O.11L1x, + 0. 102%, — 29.58%, 一 5.15.x5 
一 0.026x3 — 0.62x7 一 0.424x8 一 0.514 x3x7 . . (29) 
IRB ”外 推 63 年 
y = 52.773 — 82. 0322 一 9.513xg 一 0.776x + 35.10323 + 0.23277, (30) 


RNa Re. 外 推 63 年 F=15 
= 3.31 十 4.715x, 一 0.04x6 一 20.624x 一 1.625xs 一 0.0795x? 一 1.112x3 























十 0.00002x3 + 0.39x? 十 0.0006xsxs + 0.00108x6 (31) 
式 中 x 一 一 湖水 位 ，x: 一 一 夏 蜂 级 数 ，xs 一 一 8 月 降水 日 ，x4 一 一 7 月 降水 日 , x; 一 一 8 月 
降水 量 ， Xo 7 月 降水 量 ， Xy 7 月 温度 ， Xg 8 月 温度 。 





年 发 生 级 数 ”外 推 63 年 下 一 15 .| 
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y 一 6.048 一 0.007 xı 一 0.01x9 一 62.9 xu + 31.917 x, 一 15.347 xg 
一 0.0078x3 — 0.000008x2 + 0.000005x3 + 1.14323, 一 0.602z%, 
+ 0.322x1s + 0.0236x2%4 (32) 


























AF x: 一 一 湖水 位 ,x 一 一 上 年 9 月 降水 量 ,x 一 一 7 月 降水 日 ,xs 一 一 7 月 降水 量 ,x 一 一 6 
REKE, ru 8 HRE, rn 7 月 温度 ，xs 一 一 6 月 温度 。 
上 式 各 世代 对 1963 年 外 推 值 见 表 10。 
表 10 STR SMELL 
数值 m zl 
X a z it tt TITT 
z a EDE] ama |y y| ee ||. 5) ee ema y | ee | a 


预测 | 预测 


预测 | 预测 


预测 | 预测 预测 | 预测 


0 0 0 | 1 5.6 1.96 





从 上 述 直 线 明 归 与 曲线 角 归 在 外 推 方 面 的 结果 ,可 以 看 到 曲线 角 归 较 直 线 角 归 为 佳 ， 
则 可 能 由 于 部 分 自 变 量 与 因 变 量 之 间 的 关系 ,在 一 定 范围 内 系 呈 直线 关系 ,而 超出 此 范围 
则 有 旦 曲 线 相关 ,因此 曲 角 归 在 这 些 特殊 情况 下 郎 可 以 表现 出 让 进 归 所 不 能 表现 的 关系 。 


=. # fa 


本 区 共用 了 三 种 预 调 方法 , 郎 根据 历史 变动 规律 进行 预测 ,利用 随机 序列 与 周期 方程 
进行 预测 及 多 因 索 过 滤 授 归 预 测 法 。 第 一 个 万 法 的 依据 系 利 用 种 履 发 生 型 的 理论 ， 就 所 
得 预测 值 而 言 有 一 定 的 实用 意义 ,但 由 于 此 种 方法 在 资料 分 析 时 只 作 一 般 的 统计 处 理 , 主 

要 系 丛 观 分 析 , 因此 , 难 葛 受 主观 成 分 影响 ,尤其 当 预 测 值 发 生 在 举 间 的 “拐点 "时 ,或 对 特 
大 或 极 小 发 生年 进行 预测 时 , 可 能 有 较 多 的 误差 。 

第 二 个 方法 应 用 了 时 间 序 列 概 念 ,考虑 到 长 短 周期 特性 与 环境 因素 的 随机 影响 ,因此 
分 析 程 度 比 第 一 法 深入 ,在 开始 试用 此 项 方法 时 ,人 鲁 就 下 列 几 种 假设 进行 比较 1) 分离 时 
只 包含 趋势 项 和 周期 项 ，2) 只 包括 随机 项 ，3) 包 含 三 项 ， 但 周期 项 在 随机 项 之 前 ，4) 项 
数 与 3) 同 ,但 随机 项 在 周期 项 之 前 ,所 得 结果 以 4) 最 好 , 理由 可 能 有 二 点 ,第 一 在 飞 蝗 发 
生动 态 中 ,周期 性 不 十 分 明显 ,第 二 实地 和 经 验 说 明和 名 年 间 的 随机 影响 较 大 ， 将 随机 项 提 在 
周期 项 之 前 ,可 雹 大 其 作用 。 

另外 ,在 谐 波 分 析 结 果 中 , 5 个 周期 比较 明显 ， 除 25 年 长 周期 外 ,其 它 4 个 周期 约 小 
于 12 年 ,此 恰好 给 予 周 期 方程 中 假设 12 年 之 内 的 周期 以 理论 说 明 。 

Bye = f@) + ED) + g(t) 方程 中 分 离 因素 的 排 烈 次 序 ， 亦 可 看 出 影响 东亚 飞 蝗 发 
生 的 因素 ,随机 影响 与 周期 影响 具有 同样 重要 的 意义 , 偏 一 则 均 可 影响 预测 值 与 内 播 误差 
的 效果 。 叉 在 蝗虫 等 级 预测 中 使 用 平稳 随机 过 程 的 理论 谷 系 一 KBR, 三 稳 随机 过 程 是 
具 前 发 展 比较 迅速 的 数学 分 支 之 一 ,各 方面 已 得 到 广泛 应 用 。 

多 元 直线 迎 妇 方程 是 当前 最 常见 的 一 种 预测 式 ， 此 法 的 缺点 系 各 因素 与 种 从 的 关系 
都 是 简 相 关 , 而 名 略 变 量 本 身 相 互 的 影响 ,以 及 自 变 量 对 因 变 量 的 综合 影响 ， 多 因素 深 步 
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肖 归 法 避 竟 了 此 点 ,在 给 定 信 和 度 下 ， 权 衡 利 狗 或 轻重 抛弃 若干 变量 。 若 与 第 一 ,二 两 种 方 
法 比较 ,前 二 者 均 系 由 历史 资料 直接 进行 外 推 , 未 郑 虑 任何 具体 变量 的 作用 ， 后 者 则 考虑 
了 可 能 在 某 一 具体 时 间 发 生 的 某 些 因素 的 作用 。 

在 检验 第 三 个 方法 中 的 10 个 迎 归 式 后 ,发 现 有 9 个 通 妇 式 中 留 有 湖水 位 (变量 )， 由 
此 证 明 水 位 变化 是 影响 针 区 飞 蝗 发 生 具 有 代表 性 的 变量 。 此 外 ， 还 有 不 少 介 归 式 中 发 现 
降水 日 对 飞 蝗 发 生 有 正 的 影响 , 初 看 似乎 了 矛盾 ,因为 降水 日 增多 ,降水 量 印 可 能 增加 ,而 降 
水 量 与 飞 蝗 发 生 量 系 负 相关 ,但 若 注意 迎 归 式 中 都 出 现 有 潮水 位 或 降水 量 的 负 影 响 , 则 可 
以 理解 在 降水 量 一 定 的 条 件 下 ,分 散 于 多 日 的 降水 显 较 集 中 降水 对 飞 凰 发 生 更 为 有 利 。 

在 多 因素 逐步 凶 归 预测 式 中 ,都 设 有 二 个 信和 度 的 预测 式 (F 一 1.5，F 一 2.5), 从 数学 
理论 上 考虑 , F = 2.5 信 度 应 比 F = 1.5 可 化 , 但 实际 并 不 如 此 , 因 信 度 提高 后 ,可 能 舍 去 
若干 有 关 的 变量 ,从 而 降低 精确 性 ,作者 考虑 用 二 个 信 度 获得 的 表 妇 式 ， 至 少 可 以 从 下 列 
二 证 面 进行 比较 ;第 一 比较 二 者 的 内 插 误 差 ,我 们 知道 当 信 度 提高 时 ,一 般 变量 数 减 少 ,内 
播 误差 增加 , 若 在 下 一 2.5 下 的 内 插 误 益 远 超过 F = 1.5 的 内 插 误 益 时 , 朗 可 知道 在 提高 
信 度 水 平 的 过 程 中 , 含 去 了 有 用 的 变量 ,如 此 ,一 般 朗 认为 F = 1.5 的 通 归 式 较 好 ， 反 之 ， 
F=25 下 的 结果 则 更 可 党。 第 二 依据 实际 多 验 以 及 从 飞 蜂 生物 学 角度 进行 分 析 ， 若 角 
Ie THERES, WRIA RRM EMER, GRR F = 1.5 SQA 
式 。 

在 进行 数值 预测 中 , 1960 年 后 出 现 内 插 谋 差 负 值 较 多 , 朗 预 测 值 (预报 偏 高 ， 此 可 
兹 由 于 在 数量 对 测 的 资料 中 没有 全 部 排除 防治 的 作用 ， 因 而 所 得 预测 值 可 能 相当 于 未 防 
治 情况 下 的 发 生 数 值 ,他 如 级 数 预测 中 ,由 于 已 排除 防治 的 作用 ， 故 预测 值 比 较 吻 合 于 实 
际 发 生 数 值 。 . 

又 如 ,在 试 作 1963 年 年 发 生 级 数 的 预测 时 ,无 论 圳 信 水 平 是 Ff 二 1.5 或 F 二 2.5, 所 
预测 的 结果 与 实际 发 生 情 况 相 差 误差 在 4 级 左右 ,而 用 曲线 通 归 预测 的 结果 则 大 臻 吻合， 
分 析 其 原因 系 1963 年 淮河 中 下 游 地 区 的 降水 日 突破 所 用 历史 资料 的 最 高 记录 (70 年 一 
过 ), 由 此 可 说 明 在 一 般 情况 下 ， 蘑 些 自 变 景 与 因 变量 成 看 明 归 ， 营 出 现 大 幅度 起 越 极 限 
时 , 则 宜 用 曲 通 归 ,基于 此 , 曲 迎 归 对 东亚 飞 蝗 的 发 生 数值 预测 , 亦 有 其 实用 的 意义 。 

由 于 防治 及 耕作 等 人 类 活动 的 程度 对 飞 蝗 发 生 级 数 超重 要 的 作用 ， 因 此 在 预测 中 仅 
作 一 般 性 处 理 ,可 能 出 现 预测 值 偏 大 的 结果 ,在 短 中 期 预测 中 由 于 可 用 防治 后 的 残 凰 级 数 
作为 计算 基数 , 避 喝 一 般 性 处 理 的 不 足 , 长 期 预测 中 所 依据 的 变量 值 大 多 系 预 报 值 ， 所 以 
就 有 必要 在 长 期 预测 式 中 ,根据 当时 已 知 变量 , 列 入 一 校正 系数 。 
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STUDY ON LONG-TERM PREDICTION OF LOCUST POPULATION 
FLUCTUATIONS 


Ma SHIH-CHUN, TING YEN-CHIN & Lr DEN-MO 


_ (Institute of Zoology, Academia Sinica) 


The present paper proposes some aspects and procedures for establishing the long- 
term predictive equation of locust population fluctuations. Three kinds of predictive po- 
pulation models are introduced and discussed with respect to their practical value. A 
series of locust population data of 50 years, 1913—1962, taken from Hung-Tze Lake, a 
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typical locust region, has been used to analyse the regularity of locust population dyna- 
mics, from which the effective factors have been selected to build up the predictive 


models. 
1. Analysis of Outbreak Regularity 
1) The running equation, E = — = 2npq _ is employed to test the indepen- 


Vapa — 3pq) 
dence on outbreaks of one-another. It is found that some connections are present be- 
tween two neighbouring outbreaks. This means that each outbreak would be considered 
as a continuous part of the population fluctuation sequence which is influenced by the 
resultant effect of both locust bionomics and environmental factors. 

2) Results obtained about the influence of climatic factors, by the analytical. 
methods of contingence table and conditional probability, show that a dry climatic year is 
more suitable to the locust population growth than wet years, and a dry summer is par- 
ticularly favorable as compared with spring. Such climatic condition may be expressed by 
the climatic index shown as follows: Precipitation saturation deficit < 3 and 
<y = 0.95e°!22, (y---Mean monthly temperature during the locust abundant season) 

Besides climate, some other environmental factors have also been found to be im- 
portant in relation to locust population fluctuation or more useful for establishing the 
predictive equations, such as the monthly mean temperature of January, lake water level: 
and the mortality of overwintering eggs for the summer generation, the flooding area of 
lower lands and lake water level for the fall generation. 

3) In analyses made on periodicity of locust population dynamics, three periodic 
intervals are evident, and these are 2—3 years, 5—6 years and approximately 25 years. 
In addition, with the aid of the analysis of integral curve method, a periodicity of 25- 
year has also been revealed, and it is proved in accordance with the periodicity of dry 
months from May to August. 

4) The influence of crowding is another biological factor which has to be con- 
sidered in locust population fluctuation as exhibited by diminution of ovariole number and 
the amount of eggs. The interrelationship between them has been calculated based on a 
series experiments by using the transcendental equation which gives a well fit. 


2. Prediction equations 


The predictive models are established by the following three methods: 

1. By direct extrapolation method on the basis of 

1) the. analysis of population dynamics of the fluctuation form: First is to analyse 
the 50-year data by using .the method of transition probability matrix of Markov Chain, 
then it gives two mattices in which we may find the higher transition probabilities oc- 
curing among the series of population abundances as follows: 2 —> 3, 4 — 3, or 4 一 7， 
-> 8. 

2) the characteristics of population fluctuation sequence: It is indicated that if the 
value of population abundance deviating from the mean level (y = 4.141) is plus 
or minus 2, the fluctuation tendency will reverse, for example, when x; < xi < Xir 
then x43 < x; when x; > xi > xip then xj43 > Xita, generally xj43 < 4.1. 
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2. By the methods of random sequence and periodical equations based on environ- 
mental factors and population bionomics: the following model is used 


y, = f) + EG) + gce) 


Cyu) **. predicted value of population abundance, f(t)-++ tendency factor, E(t) --- 
random factor, g (1)..:. periodical factor) 
The estimation of these three factors dissolved are 


f*(2) = 4.475 一 0.00842 + 0.38 sin (Z e) 
ê* (e) = —0.1186E(# — 1) 一 0.14526( — 2) 
g*(e) = g(t) 十 ga(t) +eee + gnt) 


3. By stepwise multiple regression method: Equations are established based on 
numerous variates, many of which have been sifted under given confidence levels (F=1.5, 
2.5) before they are taken in the predictive models. 

1) Linear regression equation 

y = by + bixi + bira + tte + bax, 
y = ba + dix, CF test) 
y= bg + bix, + bix CF test) 


For population abundance (y) of year 1963 
y = — 26.89 — 0.742 x2 十 0.326x6 十 0.56441 + 0.859% 


2) Curvilinear regression equation 


ll 


y = ay cb ayy + axt H tte + aaa t bye + bait e 
see + b, 十 Cie yXo + C13%1%3 t- + C yjX Xj 
y = ag + aye, + bix? (F test) 


ag + aye, + apt, + bix? + by? + cori CF test) 


For the year of 1963 


y = 6.048 — 0.007 x1 一 0.01xo 一 62.9X1 十 31.917x1» 一 15.347x13 
一 0.0078 好 一 0.000008x} 十 0.000005x2 + 1.143x3, 
一 0.602x3 十 0.322x33 十 0.0236x2%4, 


It is evident that the above predictive equations based on three different methods 
represent some characteristics of the locust population fluctuation. Each of them can be 
applied independently, but the median value obtained after checking with one another has 
been proved on practice, to be more accurate. 

The term of “population abundance” used here refers to the summation of atea and 
density of locust during abundant season. l 


